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Abstract 

The weak coupling instabilities of a two dimensional Fermi system are investigated for the 
case of a square lattice using a Wilson renormalization group scheme to one loop order. We focus 
on a situation where the Fermi surface passes through two saddle points of the single particle 
dispersion. In the case of perfect nesting, the dominant instability is a spin density wave but 
d-wave superconductivity as well as charge or spin flux phases are also obtained in certain 
regions in the space of coupling parameters. The low energy regime in the vicinity of these 
instabilities can be studied analytically. Although saddle points play a major role (through 
their large contribution to the single particle density of states), the presence of low energy 
excitations along the Fermi surface rather than at isolated points is crucial and leads to an 
asymptotic decoupling of the various instabilities. This suggests a more mean-field like picture 
of these instabilities, than the one recently established by numerical studies using discretized 
Fermi surfaces. 

PACS: 71.10.Fd, 71.10.Hf, 74.72.-h 

1 Introduction 

Most of the unusual properties of the superconducting cuprates are likely to be linked to the 
quasi-two-dimensional nature of their electronic structure close to the Fermi energy. Therefore 
certain (single-band) two-dimensional (2D) models of interacting electrons may be able, in 
principle, to account for at least part of the anomalies observed in these compounds jl|, |^ . 
Unfortunately, even very simple models, such as the 2D Hubbard or the 2D t— J model, have so 
far resisted a rigorous analysis. Moreover, the available numerical studies are not yet conclusive 
enough for making definite predictions for, e.g., the zero-temperature phase diagram of these 
many-electron systems. 



One of the major difficulties is that in the cuprates the bare couphngs between electrons, 
for instance the parameter U of the Hubbard model, are large, i.e. of the order of the band- 
width. Therefore it is not clear whether a ground state consisting of occupied Bloch orbitals 
with energies below ep is a good starting point or whether one has rather to think in terms 
of configurations of singly occupied and empty sites (doped Mott insulator). Actually, the 
successful analysis of the insulating phase in terms of the Heisenberg model suggests that the 
Mott insulator is the appropriate reference state Another difficulty is that fluctuations 
(both thermal and quantum) are strong in two dimensions so that mean-field approximations 
cannot be trusted. 

In this paper we deliberately choose the limit of weak bare couplings, keeping in mind 
that this parameter range may miss completely some important characteristic aspects of the 
region of strong bare interactions. Nevertheless, it cannot be excluded that certain properties 
are qualitatively the same over the whole range of (bare) couplings, as is the case for the ID 
Hubbard model, a Luttinger liquid for all positive values of U and all densities except n — 1 



The most clear picture of two-dimensional interacting electrons has been obtained for the 
2D jellium model with its circular Fermi surface, using a Wilsonian renormalization group 
(RG) approach |0, ^ . A series of rigorous studies has shown that the Landau Fermi liquid 
theory is stable at not too low temperatures [^, i.e., above the critical temperature for 
Kohn-Luttinger superconductivity [|12|. Other instabilities do not occur. 

Electrons hopping between the sites of a square lattice yield a spectrum that differs in 
two respects from the parabolic spectrum of the jellium model. First, the spectrum exhibits 
extrema and saddle points in the Brillouin zone. General considerations imply that there are 
at least two saddle points and two extrema (one maximum and one minimum) . Obvious points 
are Po = (0, 0) and Q = (tt, tt) for the extrema and Pi = (vr, 0) and P2 = (0, n) for the two 
saddle points, but more complicated patterns are also possible. The density of states has a 
logarithmic van Hove singularity at the saddle points, in strong contrast to the constant density 
of states of the parabolic spectrum of the jellium model. The second difference is the curvature 
of the lines of constant energy. These are circles in the case of the jellium model, whereas in 
the case of the square lattice one can easily find portions with almost vanishing curvature. In 
fact, for the tight-binding model (with hopping restricted to nearest neighbor sites) the Fermi 
surface for the half-filled band is a perfect square. 

RG calculations for a model where the Fermi surface contains fiat portions have been per- 
formed by various authors |l3| , [l^ , [l^ . They agree in that a d-wave superconducting insta- 
bility occurs for repulsive interactions, due to the coupling of particle-particle and particle-hole 
correlations. 

Our main emphasis is on the effect of van Hove singularities. We will consider in particular 
the case where the Fermi surface passes through saddle points ("van Hove filling"). Early 
scaling approaches to this problem ^ focussed on the interactions between electrons 

at the saddle points, by treating these points in analogy to the two Fermi points of the one- 
dimensional electron gas In this work we show that, indeed, the logarithmically dominant 
RG flow at low energies is controlled by the neighborhood of the van Hove points. However, 
in contrast to the one-dimensional case where the scattering processes can be characterized 
in terms of a few coupling constants connecting the two Fermi points, in two dimensions 
the effective couplings are functtons of incoming and outgoing momenta, even if these are 
restricted to the Fermi surface. We find that this functional dependence plays a crucial role in 
the asymptotic decoupling of competing instabilities. A step in this direction has already been 
made in the parquet approach of Ref . |2l[ . 



When the Fermi level is at a van Hove singularity the system is not renormalizable in the 
traditional sense of field theory. Nevertheless, electrons near a van Hove singularity have been 
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treated by applying the field theory formalism to the particle-hole sector [g2| . No mixture with 
particle-particle diagrams can be treated within this formalism. The Wilsonian RG used here 
does not assume renormalizability and may be applied without constraints. 

A numerical scheme for calculating the complete flow from the bare action of an arbitrary 
microscopic model to the low-energy effective action as a function of a continuously decreasing 
energy cutoff A has been presented by Zanchi and Schulz |23[ . Unfortunately, in order to carry 
out the RG calculations it appears to be necessary to resort to a number of approximations, 
which are justified only at the final stage of the RG fiow, where A is much smaller than the 
bandwidth. Nevertheless, the application of this method to the Hubbard model near half filling 
does provide an appealing picture, namely a transition from an antiferromagnetically ordered 
ground state at half filling to a d-wave superconductor upon doping This result has 

been confirmed by Halboth and Metzner using a similar approach [M. Recently, numerical 
RG calculations have brought up two additional phases, one with a deformed Fermi surface 
(Pomeranchuk instability) 123] and one with suppressed uniform spin and charge susceptibilities 
("insulating spin liquid") p6[ . In all these calculations the proximity of van Hove singularities 
plays an important role, together with approximate nesting. 

Our analytical approach is complementary to these numerical RG calculations. We start 
from the same equations and analyze the flow in the limit of small A. We focus on the system 
at the van Hove flUing, where we take only the leading order in A into account. The asymptotic 
regime of small A can only be reached, if the initial coupling is sufficiently small. In this sense 
our approach is limited as compared to the numerical studies. On the other hand, the numerical 
methods suffer from the need to replace the continuous Fermi surface by a discrete set of points. 

In our approach - as well as in previous RG calculations carried out to one loop order - 
self-energy effects are neglected. While this can be easily justified for the jellium model, the 
argument is more subtle in the case of lattice fermions. In fact, the second-order contribution to 
the self-energy is infrared divergent in the case of the half-filled nearest-neighbor tight-binding 
band. Nevertheless, we find that also in this case self-energy effects are of subleading order in 
A, provided that an instability (superconducting or density wave) occurs. 

The analysis of the dominant parts of the RG equations is sufficient for establishing a rich 
phase diagram for the nearest-neighbor tight-binding band with a nested Fermi surface, whereas 
in the non-nested case subleading contributions are crucial. We also point out the difficulties 
of including consistently those subleading terms. 

The paper is organized as follows. In Section ^ we define the effective interaction for a 
general model of interacting electrons in two dimensions and show how to relate it to correlation 
functions. An exact RG equation and its one loop approximation is presented in Section^ both 
for the effective coupling function and for generalized susceptibilities. The close connection 
between this approximation and parquet diagrams is also briefly discussed. In Section ^, we 
analyze the one loop equations in the limit of small energies. First, we review the situation 
of a parabolic electron dispersion with a circular Fermi surface where we recover the standard 
result of a dominant flow in the BCS channel. Then a Fermi surface is considered which passes 
through van Hove points without being nested. For attractive interaction superconductivity 
again dominates, whereas it appears to be difficult to keep track consistently of all leading 
order terms for repulsive interaction. The case of a half filled nearest-neighbor tight-binding 
band is discussed in the remaining Sections In Section ^ the renormalized couplings are 
classified according to both the location of momenta with respect to the van Hove points and 
the channels characterizing the different instabilities. It is argued that to leading order there is 
no mixing between superconducting, charge and spin instabilities except from momenta very 
close to the van Hove points. A simple way of disentangling this special behavior at the van 
Hove points and the generic behavior elsewhere is presented in Section ^ and contrasted to an 
earlier approach where the momentum dependence was altogether neglected. The asymptotic 
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behavior of the RG flow aUows to draw a phase diagram including superconductivity, density 
waves and flux phases, depending on the values of the bare couplings. This phase diagram 
agrees with symmetry considerations linking the various order parameters, as shown in Section 
1^ A brief summary is presented in Section ^. 

2 Effective interaction 

We consider a system of interacting electrons on a two-dimensional square lattice. The single 
particle states are labeled by a lattice momentum k, which is defined only modulo a reciprocal 
lattice vector, and a spin index a =tii- In the functional integral formalism the system is 
described in terms of Grassmann fields il^ak, where k = (fco,k) is a 2 + 1-dimensional variable, 
which includes the Matsubara frequency fco [^. The Fourier transform is defined as 

V.fc = mr^^^ rdre^'=°^^e-"--V.r(r), 
Jo 

where /3 is the inverse temperature and V the volume of the system. In the calculations we will 
take the limit P,V oo. The action is of the form 

S[V'] = ^^.fc(«fco-Ck)V<.fc-VF[V], (1) 

where = Sk — M is the single particle energy relative to the chemical potential. The free 
electron propagator is 

C{k) = -r^- (2) 

W can be a general short ranged two-body interaction, with a coupling function g(k\,k2,k-i,ki), 

Xg(fcl, k2, k3,k4) ^ ^ak3a'k2'^cy'k3^<Tki- (3) 
a, a' 

Note that although we write formally 5k;^+k2,k^+ki, the momenta are only conserved modulo a 
reciprocal lattice vector. 

The function g should satisfy all the point symmetries of the square lattice. In addition, 
we require permutation symmetry g{ki, k2, ka, k4) = g{k2,k\,ki,k3). Finally, from time re- 
versal symmetry and the behavior under complex conjugation one finds giki, ^2, fea, fe4) = 
g(fc4, fcs, fc2, fcl) = g{ki, k2, fca, ^4), where fc = (—fco, k) and g is the complex conjugate of g. 

There is no symmetry with respect to the operation X(;(fci, fc2, fca, fc4) :— g{k\,k2,ki,kz) 
which exchanges only two of its arguments, but g can be separated into a symmetric part 
g^ = i(l -I- X)g and an antisymmetric part y'^ = i(l — X)g. These couplings describe 
scattering of singlet and triplet pairs, as becomes clear if we write the interaction as 

W[lp] = ^ 5fci+fc2,fc3+fc4 I 5^(fcl,fc2,fc3,fc4)0s(fc2,fcl)<^s(fc3,fc4) 

+S(^(fel, fc2, fca, fc4) ^ 0Q(fc2, fcl)<^a(fc3, fc4) >, (4) 

Q=0,±1 
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with 

(^s(fc,fc') = 75 (V'fctV'fc'i - ^fci^fe'r) • (5) 
The system is completely described by the partition function with source term 



where we used the short-hand notation (x, tp) := X/o-fc Xcrk^ak and the normalized Gaussian 
measure is defined by 

In particular, all connected correlation functions are obtained as functional derivatives 

52'MogZM 



(V'l • • • V'nV'n+l • • • V'2n)c = 



*72n ■ ■ ■ 5qn+i5qn ■ ■ ■ Srji 



(8) 

r]=0 



where we have written tpi instead of i/'fcjo-;- 

The propagator is singular on the manifold kg = — called the Fermi surface. The 
infinities encountered in a naive perturbative treatment of the action (|^) are avoided in the RG 
procedure which endows the bare propagator with an infrared cutoff A, 

CA(fc) = e(|Ck| -A)C(fc), (9) 

where is the Heavyside step function. The quantity that separates high- from low-energy 
degrees of freedom is therefore Although a more canonical choice would be y/kg + ^■^(p), 
we have chosen the frequency-independent cutoff in order to simplify the calculations. 
One can now define the effective interaction 

WaM = - log /dMcJ^]e-'^'*+'^' (10) 



which depends on a Grassmann field x- Note that the integration with respect to dfic/^ [tp] is 
perfectly defined, although C^^ is not. This can be seen most easily in the expansion of VVa[x] 
in terms of Feynman diagrams. The evaluation of these diagrams involves only Ca and never 
Whenever C^' appears in an intermediate step of a calculation (see below), it may be 
regularized by replacing the zero in the Heavyside function by an infinitesimal number. 

Wa has a twofold interpretation. On the one hand, we can restrict the field to the low 
energy degrees of freedom ■!/>< ka = 0(A — \^k\)tpkrT- The object 

Sf[^<] = i^<,C-'i,<)-WA[^<] (11) 

corresponds then to Wilson's effective action, which describes the system in terms of tp^ only. 

On the other hand, Wa is the generating functional of amputated connected correlation 
functions with infrared cutoff A because of the identity 

log^'AW = -iv,CAV) - Wa[Cav], (12) 
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where Za is given by with C replaced by Ca- In particular, the quadratic part of Wa is 
related to the self-energy Ea by 



SXakSXak 



x=o 



EA(fc) 



l-CA(fc)EA(fc)' 



(13) 
(14) 



where we have used the following identity for the full electron propagator G^ik) = —{"ibaki'crk) a = 
CA(fc)(l-CA(fc)SA(fe))"V Therefore in the case |^k| < A the right-hand side of Eq. ([l4) simply 
becomes EA(fc). 

Similarly the quartic part of Wa is related to the one particle irreducible vertex Fa, defined 



by {^<.fciV'<.'fc,V'a'fc3V'<Tfc4)c,A = {pvy^vT (fci 



Eq. (0 



we find 



<5'Wa[x] 



5XcTki5xa'k:i5xa>k2^<ykt 



x=o 



, ^4) Y\i-^GA{ki). In fact, differentiating 

4 

'{Tp^kiTpa'k2'^<y'k3tp,Tki)c,A C]^^ikt) (15) 

1=1 

,fc4) 



rr (fci, 



(16) 



The quartic part of Wa is of the same form as 
gA{k\, . . . ,ki). For |^ki | < A we find therefore 



/3^nr=i[i-t^A(fcosA(fcO] 

with an effective coupling function 



,k4) = -Skj^+k2,k3+ki{'i-- S„„,X)gA{ki, . . . ,k4). 



(17) 



qa is equal to a connected amputated correlation function if all \£,k-\ > A and a one particle 
irreducible vertex in the opposite case, qa is therefore not continuous at \£,k-\ = A. A formal 
and non-perturbative proof of these relations was given by Morris for a bosonic field theory 
but the generalization to fermions is straightforward. He has also shown that Ea and Fa are 
continuous at \ — A, in contrast to (?a. 



3 RG equations 

3.1 RG flow of the effective interaction 

The effective interaction satisfies the following exact RG equation |l^, ^ 

dA ^ dA ' ' SxakSxak ^ dA ' ' Sx„k Sx^k 



(T,k 



The derivative in ^Ca restricts the propagator to an infinitesimal energy shell (A, A -I- dA) 



and Eq. dig ) therefore describes the effect of integrating out the modes of that energy shell. 
This equation was first derived by Polchinski in the context of a scalar field theory [po| . 

The strategy is now to start with the bare interaction at a cutoff Ao and to use some 
truncation of Eq. ( [l8| ) in order to compute approximately the effective action in the limit 
A — > 0. We will follow here Zanchi and Schulz 123] who proposed to develop Wa up to order 
six in the fermionic variables and to neglect terms of higher order. 

We will only follow the flow of the effective coupling function and neglect self energy cor- 
rections (i.e. the quadratic part of Wa). Self energy corrections to the single particle Green's 
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PH2: 





ks 
k2 




Figure 1: Diagrams contributing to the renormalization of the coupling function gA(fci, . . . , A;4 



function have three main effects. First they change the shape and location of the Fermi sur- 
face, second they modify the properties of the single particle dispersion (the Fermi velocity) 
and third they lead to a reduction of the quasi particle weight. It has to be chec ked from case 



to case whether such corrections can be safely neglected or not (see Sections 4.1 



4.2 and 



Terms of order six are not present in the original interaction but they are produced by the 
RG flow. Their effect is then to renormalize the effective coupling function qa- The result is a 
closed one-loop equation for the coupling function g^iki, . . . , k^) |P3[ where j^k; ] < A. It reads 



dA 



,ki) = PP-HPH1 + PH2 



PP 



PHI = 



-y 



d[CA(p)CA(go)] 
dA 



g-Ko {ki,k2,qo, p)gAg (p, go , ^3 , ^4 



1 y^ d[CA(p)CA(gi)] 
PV ^ dA 



gAiiki,qi,k3,p)g^^{p,k2,qi,k4), 



(19) 



= ''^'''^dA'^''^^ [-2g^^(fc.,p,g.,fc4)gA.fe.fe,fc3.p) 

p 

+5a2 {ki,P,k4,, q2)gx^ (92 , ^2 , fcs , p) -f g-^_^ (ki,p,q2, ki)g-^_^ (q2 , , p, ^3 )] , 

where qo = ki -\- k2 — p, qi = p + k^ — ki, q2 = p ki — k2 and K = Max{|5p|, |^qj}. 

Note that unlike Eq. (^8|) which is valid for any choice of cutoff function Ca, this truncated 
equation assumes a sharp cutoff as in Eq. (^. Its representation in terms of Feynman diagrams 
is shown in Fig. ^ 

One of the two internal lines stands for a propagator Ca and the other for its derivative 
^Ca. This amounts to a loop integration over an energy shell A < |^p| < A + dA with the 
restriction that the second propagator is in the high energy regime |^qj > A. The wavy lines 
in the diagrams represent effective couplings of an earlier stage of the RG flow, when the 
energy shell A was integrated out. 
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Eq. ( |l9[ ) is not local in the variable A. Since this is not very convenient, it was proposed 
to develop Eq. into Wick ordered polynomials of the fermionic variables instead of 

monomials as it was done above. Wick ordering with respect to the low energy propagator 
Da = C — Ca results in the same one-loop equation as above but now all the couplings 
are evaluated at the actual RG variable A and — d [CA(p)CA(g)] /dA has to be replaced by 
d [DA{p)DA{q)] /dA, i.e., the energy of the second propagator is now restricted to be smaller 
than A. 

Recently a third scheme was established where the equation is local in A and the 
second propagator is still in the high energy sector. It is obtained by considering the generating 
functional of the one particle irreducible vertex functions rather than the effective interaction 
defined by Eq. ([lo|), an idea originally due to Wetterich |^l|. 



3.2 Order parameters and generalized susceptibilities 

The RG formalism can be used to calculate the linear response to an external field p^ . 
We focus on the prominent instabilities of the half-filled nearest-neighbor tight-binding band 
by adding the term 

_ 1-0 _ 

W'[h,'tP]^~h dTO{T)~h drO{r), (20) 
Jo Jo 

to the action, where h is the external field /i G C and O is a characteristic order parameter. 
As a first example we consider superconductivity, described by the order parameter O = 
/(k) ipikipi-k with a form factor /(k) depending on the symmetry chosen (s-wave, d-wave, 
etc.). The effective interaction, defined by Eq. ( p^ ) with W replaced by -I- W' , is now 

Wa[/i,V]= WaWI^^o- ^^Xl^'^('')'^TfcV'i-fc+h.c.^ -XAPV\hf (21) 

plus terms of higher order in h and i/'- The coefficients RA{k) and xa are actually defined by 
Eq. (^l|). Their initial values are RAo{k) — /(k) and xaq = 0. The susceptibility xa gives the 
linear response with respect to the perturbation (^^ of the system with an infrared cutoff A. 
Using Eq. (|2l[) and the definition (hol) we get 



1 rwA[/i,V] 

Xa 



PV dhdh 



1 d{0)A 



(22) 



Xa is therefore interpreted as a susceptibility at temperature T = A [^. In fact, in finite 
temperature perturbation theory the temperature acts like an infrared cutoff. 

One can now derive one loop RG equations for i?A(fc) and xa following the same procedure 
as for the effective coupling function qa- We find 

k 

^^A(fc) = -^^^[Ca(p)Ca(-p)] <?r^(fc,p)i?A(p), (23) 
p 

where the BCS coupling function Qa"^^ {k,k') :— gA{k, —k' , k') describes the scattering of 
Cooper pairs with zero total momentum. Note that the function i?A(k) has the same symmetry 
properties (s-wave, d-wave etc.) with respect to the lattice point group as the initial form factor 
/(k). 
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Figure 2: The pattern of charge (spin) currents along the bonds of the square lattice in a char 
(spin) flux phase 



At half filling, instabilities related to the nesting vector Q — (tt, tt) are also natural candi- 
dates. They are described by the order parameter 

= J2 /-(k) ^-k+Q V<Tk = e""^ /4r - r') i;.ri>.r', (24) 

k,a r.r',(j 

where / is the Fourier transform of /. We assume that f<y(k + Q) = fcr(k) so that O is 
Hermitian and ft £ R. 

The simple choices f<7(k) ~ 1 and /^(k) = a yield, respectively, charge and spin density 
waves. If fcr{k) is a function with ,,2 -symmetry and is even (odd) in the spin index, then 
the nearest-neighbor-terms in Eq. (UM yield circular charge (spin) currents flowing around the 
plaquettes of the square lattice with alternating directions (see Fig. ^ . These four charge and 
spin instabilities have been discussed a long time ago in the context of the excitonic insulator 



We call the phase with circulating charge currents the charge flux phase (CF) it is 

sometimes also called d-density wave, charge current wave or orbital antiferromagnetism. The 
charge flux phase (CF), closely related to the concept of the chiral spin liquid [K4|, still plays 
a prominent role in the strong coupling 5(7(2) theory of the t — J model [p5|. Recently it was 
proposed to be the competing order parameter to d-wave superconductivity and responsible for 
the pseudo gap phase of the cuprates |Q. 

The phase with circulating spin currents is called the spin flux phase (SF). Other names 
encountered in the literature are "spin current wave" or "spin nematic state" (because it is 
a state with broken rotational symmetry and unbroken time reversal symmetry). The low- 
temperature thermodynamics of both the charge and spin flux phases have been investigated 
within a mean held theory in Ref. ^t\ . 

The effective interaction in the presence of such a perturbation is of the form 

WA{h)[i;] = >Va(0)[V'] - 2/i^i?,A(fc)^^fc+Q Vafc -XAPVh^ + (25) 

k.a 

where Q :— (0, Q) — (0, 7r,7r). The one loop equations involve the scattering processes 
g^'ik, k') :— Xg^{k, k') :— g{k, k' + Q, k' , k + Q) with a direct or exchanged momentum transfer 
of Q. They read 



= -S^^E + I^-a(p)^ (26) 



_d 

dA'^^' ~I3V ^ dA 

p,cr 
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p 



where A = Max{|^p|, |^p+q|}. One can decouple Eq. (|2q ) into spin- and charge sectors by 
writing X = X*^ + X* 



d 



c.s 



dA i3V ^ dA 

^i?l'=(fc) = ^^A[CA(p)CA(p + Q)]4-(fc,p)i?2^(p), 

where i?'^'" := i?T ± -Ri, 5" := 2^'' - and = -p^. 

3.3 Relation to parquet diagrams 

Instead of calculating the effective interaction via the Polchinski equation (^), one could try 
a naive perturbative expansion (starting from Eq. (p^). 

5A = PAo + ^1 (A) + A2(A) sio + ■ • ■ 

The problem with this series is that the coefficients An (A) diverge as A ^ 0. Their asymptotic 
behavior is given by j4„(A) ~ i"(A), where in general /(A) = \ logA| except for van Hove filling 
where Z(A) — log'^ A. If the bare interaction gtva is small but gt^^li^K) ~ 1, it is a reasonable 
approximation to sum the whole series, treating every term to leading logarithmic order in A. 
This amounts to summing the so-called parquet diagrams (see §) for a detailed description of 
this method). 

It is interesting to note that these diagrams are generated by successively integrating Eq. 
(^) and expressing the result in terms of the bare interaction -0 We conclude that solving 
Eq. ([19I ) to leading order in A is equivalent to a summation of the parquet diagrams within 
logarithmic precision. 



4 The RG flow in the asymptotic regime 

The truncation schemes applied so far to Eq. ( p^ are all of perturbative nature and only 
justified as long as the effective couplings are weak. Nevertheless, Eq. (|l^) and also its variants 
are still too complicated to be dealt with numerically or analytically. We therefore search for 
a second small parameter, in addition to the coupling strength. 

A natural small parameter is the energy cutoff A at the final stage of the RG flow, where 
the important contributions come from momenta close to the Fermi surface. At this final stage 
we approximate the coupling functions (/A(fci, . . . , ^4) by their values on the Fermi surface, i.e., 
the 2-1-1 dimensional momenta ki in gA{ki, . . . ,^4) are replaced by suitable projections on 
the Fermi surface (^k = fco = 0). We therefore write (jiA(ki, . . . , k4) instead of Qhiki, . . . , fe4), 
because the frequencies have been set to zero. The geometrical constraint of locating all the 
four momenta ki,k2,k3 and k4 = ki -I- k2 — ks close to the Fermi surface often reduces the 

^ The structure of Eq. (|is|) introduces a constraint on the energies of internal lines. For example, if a parquet 
diagram is generated by insertion of a bubble B\ into a bubble B2, both single particle energies of B\ have to be 
bigger than those of B2- The different RG schemes ]23| p^ , all generate the whole series of parquet diagrams, but 
the constraints are different. 
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number of variables even further, so that depends on two or at most three independent 
(angular) variables instead of nine. 

This restriction to the final stage of the RG flow of the effective couplings implies that we 
cannot relate the effective theory to the original microscopic Hamiltonian. More importantly, 
we have to assume that no instability occurred in the flow before we have reached the asymptotic 
region. 

In the following we identify the leading contributions to Eq. (|l^) in the limit where these 
simpliflcations are justified, i.e., in the asymptotic regime A ^ 0. The external momenta are 
assumed to lie on the Fermi surface. 

We first observe that the contributions to Eq. (|l^) are of the form 

~PV^ dA ffA(---)ffA(--0 (28) 

p 

where k — ki + k2 in diagram PP, k — k^ — ki in diagram PHI and k = ks — k2 in diagram 
PIf2. The minus sign appears in the particle-particle (p-p) diagram PP and the plus sign in 
the particle-hole (p-h) diagrams PHI and PH2. 

If we neglect for the moment the angular dependence of two coupling functions we are left 
with an analysis of the two bubbles, 

R (\k\- ^ d[CA(p)CA(fc=Fp)] 

-«pp,ph(A,A;)-- — 2^ — . (29) 

p 

We consider the thermodynamic limit and zero temperature and therefore replace 1/pV y^,_ by 
j (27r)^+i calculations. Taking explicitly the derivative with respect to A and integrating 

over the frequency po (for ko = 0) we find 



i?,,.,,(A,k)^±2 / Ji^^dgpl-A) ^^'^-^' / .r, ^ (30) 



d^p .n, I _ , ^ e(ICq|-A)e(± gpCq 

where q = k =F p. 

Within the Wick ordered scheme |^^, the first step function in Eq. ( |3o| ) is replaced by 
0(A — Ifql), since the second propagator in Eq. ( p^ is restricted to be in the low energy part. 
We have verified that this alternative scheme would not change our final results. 



4.1 Circular Fermi surface 

In order to illustrate the approach we consider the example of very low filling, where the single- 
electron spectrum is approximately parabolic Cp = ~ 1 (^^^ energies are given in units of the 
Fermi energy and all momenta in units of kp). The energy shell consists of two circles with 
radius \/l ± A. 

In this case, self energy effects are not expected to change the leading order result. On the 
one hand, the shape of the Fermi surface is flxed by the rotational symmetry of the problem. On 
the other hand, perturbative corrections to the Fermi velocity VpEa(p) and the quasi particle 
weight z = [1 + idpg'EA{p))~^ are finite as A — > 0. 

For general values of k one finds that both bubbles i3pp,ph(A,k) are proportional to log A. 
For the special value |k| = 2 we find Bpp ~ Bph ~ A~^^^. The strongest divergence comes 
from the p-p bubble at k = 0, namely Bpp ~ A^^. Correspondingly, the dominant contribution 
in the asymptotic regime A — > comes from the diagram PP. It renormalizes the coupling 
function g^'^^ {k., k') := (?A(k, — k, — k', k') which is related to superconductivity, as was shown 
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in Section 3.2, The non-locality in A of Eq. is absent if only this diagram is taken into 

account. 

The first sub- leading contribution (~ A^^''^) concerns scattering processes with momentum 
transfer 2kF- This contribution is usually neglected. We therefore conclude, in line with the 
more standard scaling analysis used in previous works Q , that the dominant flow renormalizes 
only the BCS couplings through the particle-particle diagram. 

It is worthwhile to discuss the behavior of the bubbles for small but finite values of k. For 
the p-p bubble we find in the limit |k|, A ^ 1 



1 I 2 arctan 





k| 


A+ 


k 



Bpp{A,k) = { yj^i ' ' (31) 



If we renormalize a coupling with a small total momentum k = ki -I- k2 the p-p contribution 
to its fiow is independent of k as long as A 3> |k|. In this case we replace gA(fci, . . . , ^4) by 
g^^'^ (ki,]s.4,) even if the total momentum is not exactly zero. This replacement is no longer 
justified when A is of order |k|. The flow then depends strongly on k and cannot be controlled. 
Nevertheless, there is no danger because a few renormalization steps later, if A ^ |k|, the flow 
is suppressed and the coupling under consideration no longer contributes. 

We can use the same type of analysis for the p-h bubble with a small momentum transfer, 
which renormalizes couplings that are close to the forward- or exchange scattering g^(k,k') ~ 
XffX(k,k') := c/A(k,k',k',k). We obtain 

f if jkj < A, 

i -2^1og^T if|kl>A. 

It gives a big k-dependent contribution if k is of order A. But this flow is again suppressed if 
A is further reduced. 

The results presented above for the isotropic case should remain valid as long as the Fermi 
surface is both far away from van Hove singularities and not nested. Even the presence of 
Umklapp scattering alone does not change the result. If the Fermi surface is big enough, 
Umklapp processes open up new scattering processes, but their contribution to the flow is not 
of leading order jso) . They can nevertheless influence the flow in an important way before we 
enter the asymptotic regime j26|. But at this early state of the RG flow, approximations to 
Eq. (^9]) cannot be controlled by the small parameter A. 



4.2 Generic Fermi surface at van Hove filling 

We now consider the case where the Fermi surface passes through a van Hove singularity. 
The generic situation occurs for a band structure with saddle points at Pi = (±7r,0) and 
P2 = (0, ±7r). To be specific, we consider the dispersion relation of a generalized tight-binding 
model: ~ — 2(cosfci: -I- cosfc^) + 4t' {cos cos ky + 1). The unit of energy is given by the 
hopping amplitude between nearest neighbors and the unit of momenta is the inverse of the 
lattice constant. A finite electron hopping < t' < 1/2 between next nearest neighbors has 
been included and the chemical potential is fine-tuned such that the Fermi surface contains the 
saddle points (see Fig. |^. 

We identify again the dominant terms of order A~^ in the RG equation ( |l9| ) by analyzing 
the bubbles Bpp{A, k) and Bph{A, k). Dominant terms of this order appear only if k is close to 
or close to (vr, vr). 
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Figure 3: The Fermi surface at van Hove filling: t' = 0.3 (solid line) and t' = (dashed line). 



We first study the behavior of the bubbles at small |k|. It is instructive to focus on the 
contributions of a small patch surrounding a saddle point, say the region 

Pi = {p ; VI - 2t' - 7v\ + Vl + 2t' \py\ < 2p} . 

The parameter p is small enough so that can be replaced by its limiting quadratic form close 
to Pi. We compute the restricted bubbles Bpp (A, k) and (A, k), which are defined by Eq. 
(^o|), but restricting the summation to p £ Pi. The values of Bpp (A, k) and B^^ (A, k) depend 
sensitively on ^p^+k. Both bubbles are negligible if |Cpi+k| 2> A, but for j^Pj+k| ^ A (and 
A, k^ < p^) we get 



B^^{A,k) = 



log 



4Ap^ 



and 



S,^;(A,k) = 



(27r) Vl - 4t'2 A 



(27r) VI - 4t'2 A (A + /5fc+)(A + /9fc_) 



9{pk+ - A)!—^ + e{pk- -A)^ 



pkj, 



pk- 



(33) 



(34) 



where k± ■- \^/Y~^2t' k^ ± Vl + 2t' ky\. Note that k+ fc_ = I^Pi+kj. 

There are several remarkable differences with respect to the circular Fermi surface. First 
we observe that for very small k, such that pk± <^ A, we have P^(A,k) ~ A^^log ^j^, i.e. 
the p-p bubble has a logarithm in addition to the A~^ behavior. This anomaly, due to the 
diverging density of states at the Fermi level, concerns only the renormalization of the BCS 
couplings. 

More striking is the behavior in the regime ^ A <C pk±. The p-h bubble, zero 

for sufficiently small |k|, gives a non negligible contribution there. Thus for a given small 
momentum transfer k, the p-h contributions [diagram PHI or PH2 of Fig. |l] are dominant 
over many RG iterations, in contrast to the case of Section without van Hove singularities. 
The p-p diagram is very sensitive with respect to both the size and the direction of the total 
momentum k in this regime, i.e. over a sizeable range of energy scales. We conclude that in 
a consistent treatment of the RG equations to order A~^, the renormalization of the effective 
coupling function (;A(ki, . . . ,k4) depends on the exact values of its arguments. For example 
the abovementioned projection of the momenta onto the Fermi surface cannot be justified in 
the vicinity of a van Hove point. 

One way to reduce the complexity of the problem is to keep only the terms of order 
A~^log(l/A) and to neglect all other contributions. This leading logarithmic contribution 
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occurs only in the p-p bubble for |k| <^ y/X. As a consequence only the BCS coupling function 
g^'^'^ (k, k') is renormalized within this approximation. The flow of g^^'^ (k, k') is similar to the 
one explained in detail in the next section, but without competing charge and spin instabilities. 
For k close to (7r,7r) and finite values of t' the p-p and p-h bubbles behave like but not 
like A"^ log A. 

For repulsive interactions, where the BCS flow is towards weak coupling, the approximation 
is clearly not sufficient and a consistent treatment up to order A~^ is required. Solutions for this 
case proposed earlier have neglected the sensitive dependence of 5(ki, . . . , k4) on its arguments 

©BP- 

In contrast to the case of the circular Fermi surface, self-energy corrections are not negligible 
in the present situation. In fact, they are expected to change the shape of the Fermi surface. 

For the corrections to the single particle dispersion and the quasi particle weight, one finds 
in second order perturbation theory (see for example Eq. 7 of |}4C| ) 

3poSa(p) ~ 5i„log^A (35) 
VpEa(p) ~ sio log' A ■ VpCp 

These corrections are negligible within our approximation, as long as log^ A ~ 1. An 
attractive coupling diverges at this scale. In the case of repulsive interactions however the RG 
flow can be followed to smaller energies such that ^Aq log A ~ 1. At these energy scales, self 
energy corrections have to be taken seriously. This would require to include subleading terms, 
a task that cannot be fulfilled consistently without going beyond the one-loop approximation. 

Due to the additional logarithm in the p-p bubble the theory is not renormalizable in the 
usual sense of field theory, i. e. it is not possible to send the bare momentum cutoff to infinity 
while keeping some physical correlation functions finite (even after the introduction of counter- 
terms in the microscopic Hamiltonian and of wave-function renormalization) . Gonzalez, Guinea 
and Vozmediano j2^] proposed a field theoretical RG scheme where the coupling constants for 
forward and exchange scattering are renormalized only by the p-h diagrams. In this approach 
the p-p channel is treated separately in connection with a renormalized chemical potential. It 
is however not clear from our Wilsonian approach that the RG equations will not mix p-p and 
p-h diagrams, if all contributions ~ A^^ are taken into account. 



5 The square Fermi surface 

In this and in the following Sections we consider the nearest-neighbor tight-binding model 
(t' = 0) at half filling where the Fermi surface is a square (see Fig. |^. We first identify the 
possible scattering processes which connect four momenta on the Fermi surface and satisfy 
momentum conservation modulo a reciprocal lattice vector. To simplify notation we omit the 
subscript A in ^a. In addition to the usual forward, exchange and BCS scattering, there are the 
two kinds of processes related to direct or exchanged transfer of the nesting vector Q — {n, n): 
g'^{k,'k') := Xg'^{k,'k') = 5(k,k' -|- Q,k',k-|- Q). Furthermore if three points ki,k2,k3 are 
chosen freely on two parallel sides of the square, the resulting k4 — ki-|-k2— ka lies automatically 
on the Fermi surface, giving rise to a three-parameter family 3" (ki , k2, ks). Finally, there is 
scattering of pairs with a total momentum Q, described by g^(k, k') := g(k, Q — k, Q — k', k'). 
Some examples are shown in Fig. |^. 

This list of the possible low energy processes is complete but the classification into g^ , g"^, 
S^'^'^i g'^j 9^ J 9^ S'lid (?" is not unique. For example if k and k' belong to the same pair of 
parallel sides of the square, the two-parameter families g^ , g'^, g^'^'^ , g'', g^ and g^ belong to 
the larger three parameter family g". Furthermore, two two-parameter families intersect in a 
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Figure 4: Examples of processes from the different interaction channels 



one-parameter family as ff^'^'^(k, k + Q) = g''(k, Q - k), ^^'^'^(k, Q - k) = ^^(k, Q - k), etc. 
Finally three two-parameter families can intersect in scatterings between the two saddle points 
as 5«^«(Pi,P2) ^^''(Pi.Pa) = <7^(Pi,P2). 

Due to exact nesting, £,p+Q ~ —£,p, the non interacting model has particle-hole symmetry. 
We further assume that the interaction also respects this symmetry (see Section ^). The self 
energy then satisfies the exact relation Ea(po,p) ~ — Ea(— po,p + Q), and therefore vanishes 
on the Fermi surface (po = = 0). This means that the Fermi surface is not modified by 
self-energy effects. 

For the corrections to the single particle dispersion and the quasi particle weight, one can 
derive similar relations as in Eq. (pSl), namely 

9poEa(p) < const ■ log^ A (36) 
VpEa(p) < const • gi^log^A • Vp^p. 

These corrections are negligible within our approximation, where ^Aq log^ A ~ 1 in the spirit 



of the parquet approximation (see Section 3. J ) 
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Figure 5: Plot of the p-p bubble for k = (k, k), parallel to the square Fermi surface and A = 



We now investigate the one-loop corrections to the coupling constants by an analysis of the 
bubbles (^). As a consequence of the particle-hole symmetry, the p-p and p-h bubbles are 
related by 

Bpp(A,k) = -Bph(A,k + Q). (37) 

Both the p-p and p-h bubbles are of order A^^ whenever k = n{-K, tt) + ±1) (n £ Z, k G 
R). It follows that the three parameter coupling function is renormalized by contributions 
of order A~^ from every diagram PP, PHI and PH2 in Eq. (^^. The reason is that the 
Fermi surface consists of straight lines. Namely if k is parallel to such a line, and ^p+k 
are both small for many values of p, which gives rise to a big value of the integral in Eq. 
(|^. In Fig. I we show a plot of _Bpp(A, k) for k = (k, k). For A < k < 1, Bpp(A,k) is 
well approximated by (27r)~'^A~^ log -, whereas for k <c; A it approaches its maximal value 
i?pp(A,0) = (2^)-2A-ilogi^. 

Due to this big variety of contributions of order A ^, the flow to that order appears to be 
too complicated for analytical treatment. Numerical treatments have been presented by Zanchi 
and Schulz jisj and more recently by Halboth and Metzner [^^. A parquet solution for a flat 
Fermi surface has been given by Zheleznyak, Yakovenko and Dzyaloshinskii [ [L5[ . 

By restricting ourselves to the logarithmically dominant terms of order A~ log A we can go 
a long way using an analytical approach, as will be shown now. We consider the two-parameter 
coupling functions g^'^^ , and g''. g^'^^ gets a dominant contribution (~ i3pp(A, 0)) from 
the diagram PP, g^ from PHI and g'^ from PH2. For example we have 

-^i-g^'^'^(k,k') = PP-FPH1 + PH2, 
dA y ' ' 

where PP ~ Bpp(A, 0) ~ A"^ log(l/A), PHI ~ Bph(A, k -f- k') and PH2 ~ Bph(A, k - k'). We 
may neglect PHI and PH2 for small values of the cutoff, except when k ± k' ~ Q. Similarly, 
for generic k and k', g^(k,k!) is renormalized only by the diagram PHI and (;'*(k,k') only by 
PH2. The flow equations (IL9|) are again local in A and they decouple into three independent 
identical equations for the functions g^^^ , g^ = —g^ and g"^ = 2g'' — g^ . They read 

2AA5-(k,k') = 1^5(1Cp| - A)5*(k,p)/(p,k'), (38) 
p 
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Figure 6: The Brillouin zone with the energy shell |^p| = A. It is separated into two patches Pi 
and P2 and the remaining part. 



where o stands for BCS, s or c. This is evident for g^^^ and g" since only one diagram is 
involved in their renormalization. The PH2 contribution to g'' on the other hand consists of 
three terms (the three PH2 diagrams) involving g and p^. We see with amazement that the 
special combination g'^ — 2g'^_-g^ satishes the same closed and simple equation as g^^^ and 



g" . It was shown in Section 3.2 that g^^^ is related to pairing, g" to spin instabilities and g'^ 
to charge instabilities. 

In the integration over p in Eq. (^) a large contribution comes from a small neighborhood 
of the saddle-points, due to the diverging density of states. To identify the logarithmically 
diverging contribution to Eq. (js^) we consider patches Pi = {p ; \px ~ ^| + IPal < 2p} and 
P2 = {p ; \Px\ + \Py — 7r| < 2p} of a size p <^ n/2 around the two van Hove points and separate 

the integral Epinto EpgPj + EpgPa + 5I]pGi3.z.-Pi-P2 ' ^^^''^ ^'^^ Brillouin 
zone (see Fig. |6|). We compute the weight of the patch and of the remaining part of the 
Brillouin zone, assuming that (/"(k, p) is of the same order of magnitude for every value of p. 
Comparing the values 

^E'5(l^p|-^)-2il°s(^), (39) 

y E '5(l«pl-^) = ^i°g(;^)' (40) 
PSB.Z.-P1-P2 

we conclude that the patch contribution dominates the remaining part if A ^ p^. Under the 
hypothesis that the functions (/*(k,k') are slowly varying, it is then consistent to replace Eq. 

^g*(k,k') = -E <7nk)ff?(k'), (41) 

i = l,2 

where we have set (?*(k, p) ~ (?f (k) for p £ Pi and I =_o^ log^(^^) the new RG parameter. 
Correspondingly, setting (k') « gr* for k' G Pj in ( [ll] ) leads to 

j^g]{k) = -J2 9t{^)9h- (42) 
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We are now left with the problem of renormalizing the couplings g^j, which describe scat- 
tering of particles near the saddle points. Eq. ps] ) cannot be used here because it is not valid 
for k± k' ~ Q, i.e. if k G Pi and k' £ P2. The RG flow of the couplings gfj will be discussed 
in detail in Section |^. 

Once the evolution of the as a function of I is known, one can integrate Eqs. ( ji^ ) and 
(El. Introducing ^^(k) = 3?(k) ± (k) we get 



ff±(k) = g±i^)\i=i„ ■ exp 



(43) 



and 



fl*(k,k') = 5"(k,k') 



dl [5;(k)5;(k')+/_(k)/-(k')] 



(44) 



We will see that the RG equations of g^j as a function of I yield diverging solutions at a 
finite value I = lc- Near this critical value they behave asymptotically like 



L-l 



where the constant gfj can be determined from the RG equation and a > —1. Eqs. 
( ^ then give 



5i(k) ^ A± . gim.^i, ■ (Ic - + 0{l. ^ /) 



(45) 
and 

(46) 



and 



^(k,k') 



(k,k')| 



5."(k)5^(k')| 



Mill {0,29* + c«+2} 



(47) 

where g± :— gii ± §12 and A±, B± are positive constants. The coupling function (7*(k,k') is 
diverging if g^ or < —1/2. 

The functions (?+(k) have s-wave symmetry, i.e. they respect all the point symmetries of the 
square lattice. On the other hand we see that 3* (k) is of the d^2_j,2-wave type [g-{kx,ky) = 
g-{kx, —ky) = —g'L{ky, kx)]. Thus the diverging part of the coupling function has s or d^i_y-2- 
wave symmetry. 

In the preceding calculation we have distinguished strictly between points far from the 
saddle points and those close to the saddle points. The scale which distinguishes between "far" 
and "close" is the patch size p, which was introduced by hand. 

The behavior of the overall coupling function ^'^(k, k') near the critical point depends on 
the smallest (i.e. the most negative) value of the constants g%. In Section 3.2 we will show 
that Min{§5:} = —1- In this situation the function (7*(k, k') diverges everywhere with the 
same power (/ — Ic)^^ and g*(k, k') remains a smooth function upon renormalization even if 
k or k' (or both) approach the saddle points. This justifies a posteriori the estimation of the 
diagrams (S) by the bubbles (^^ as well as the approximations of Eqs. (^) and ( ^^ , namely 
that (/*(k, k ) approaches continuously (?f(k) as k' Pi and (/f(k) approaches gfj as k — > Pj. 
If instead we would find > Min{^5:} > —1, the scattering of particles near the van Hove 
points would diverge more rapidly than that of particles at the remaining Fermi surface. This 
would mean that different regions in the Brillouin zone behave differently. The region around 
the saddle points would become strongly interacting while the remaining Fermi surface would 
remain weakly interacting. 
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In order to calculate the susceptibilities we consider Eqs. (|2^) and (|2j). The non- locality 
in the variable A of Eq. (|2^) disappears because of exact nesting (Cp+q ~ ~?p)- In the low 
energy regime A — + 0, we assume that the frequency dependence of R is not important and 
replace R{k) by i?(k). We can then explicitly perform the frequency integral to obtain identical 
equations for the charge, spin and pairing susceptibilities 

2A^X^ = -i^5(|Ck|-A)|i^o(k)|^ 

k 

2A-^R\k) = 1^5(|Cp|-A)/(k,p)i?^(p), (48) 
p 

where o = c, s, BCS. To simplify the notation we have omitted the index A in the symbols x, 
R and g. 

We treat Eq. (^) in the same way as Eq. (^8|), i.e. we take the leading contribution from 
the patches Pi,P2 around the two saddle points and assume -R*(p) ~ Rt for p G Pi. We get 

^^Rl = -glRl, (49) 

where R% := i?* ± _Rf and x" —: x+ + X-- K 9± is diverging asymptotically like g±/{lc — I) 
and if g± < —1/2, the corresponding susceptibility diverges with a critical exponent + 1: 

xi~(/e-0'^±+'~(A-A.)'^±+^ (50) 

On the other hand g± > —1/2 leads to a finite value of the susceptibility. 

We thus naturally identify six possible instabilities. In each of the charge, spin or pairing 
sectors the form factor can be either even or odd under the exchange of the two van Hove 
points (i.e. Ri = R2 or Ri = —R2)- In the pairing sector x+'~'^ clearly corresponds to s- 
wave superconductivity (sSC) and X-^^ to ci^2_y2-wave superconductivity (dSC). The charge 
and spin density waves are related to the susceptibilities xX s-nd x+i respectively. The two 
remaining susceptibilities X-" correspond to a form factor with d^2_y2-wave symmetry in the 
charge and spin sectors. They describe the tendency towards the formation of charge or spin 
flux phases (see Section 3.2). 

These six instabilities of a system with two van Hove singularities have been discussed long 
ago by H. J. Schulz [Q. Here we have shown that they appear naturally in the very weak 
coupling limit of Wilson's renormalization group. 



6 Flow of the couplings between saddle points. 

We now return to the renormalization of coupling constants for scattering processes both within 
and between saddle point patches. We consider the nearest-neighbor tight-binding model (t' — 
0) as in Section ^. The most simple assumption is to treat the saddle points in close analogy 
to the Fermi points of a one-dimensional system where there are just four types of scattering 
processes, one restricted to the region of a single Fermi point, the other three involving both 
right and left movers (forward, backward and Umklapp scattering). This one-dimensional 
scenario with four coupling constants gi, . . . ,g4 implicitly assumes that the more detailed wave 
vector dependence in this region is irrelevant. While this appears to be true for one- dimensional 
Fermi systems, where going away from a Fermi point means leaving the Fermi surface, we will 
argue that in the present case the functional dependence of the couplings is relevant in the 
neighborhood of the saddle points. 
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6.1 "One-dimensional" solution 

In an early contribution to this subject H. J. Schuiz assumed that the couphng function 
g{\i.i, . . . , k4) takes only four different values according to how the momenta ki, . . . , k4 are 
distributed over the two patches, namely 

{gi ; ki , ks G Pi and k2 , k4 G P2 
g2 ; ki, k4 G Pi and k2, ka G P2 . 
ff3;ki,k2GPi andk3,k4GP2 ' 
34 ; ki , . . . , k4 G Pi 

Note that we have interchanged g2 and 34 as compared to Ref. |f^. The parameters gj- = 
gii ± gi2, which control the various instabilities (o = s,c, BCS), are readily expressed in terms 
of the coupling constants gi, . . . , 34, 



SDW/SF 
CDW/CF 
sSC/dSC 



9± = -92 T 93 
9± = 2^1 - 52 ± ffs 
ffl^^ =<?4±53. (52) 



gi has to be renormalized by the diagrams PH2, because the direct momentum transfer 
ks — k2 is close to the nesting vector Q, but the other contributions coming from PHI and PP 
are negligible. Similarly 512 is renormalized only by PHI and 514 by PP. The remaining coupling 
g3 on the other hand gets leading contributions from all three channels, PP, PHI and PH2. 

The RG equation is obtained by locating the external momenta ki , . . . , k4 in Eq. (|l|) 
exactly at the saddle points Pi and P2 and restricting the sum over p to the two patches 
Pi,P2. The resuh is 

^51 = 231(32-51) 

d 22 
^52 = 52 + ffa 

^53 = 233(252-51-54) 

^54 = -53 - 54- (53) 

For most initial conditions the numerical solution of these equations diverge asymptotically 
like in Eq. (^) with coefficients gi satisfying 

51 = 251(52-51) 

52 = 52 + 53 

53 = 253(232-51-54) 

54 = -53 - 54- (54) 

Eq. ds^ ) has many solutions, but the ones which are relevant for the divergences of Eq. ( js^ ) 
are 51 = 0, 52 = —54 = 1/6 and 53 — ±\/5/6, depending on whether the initial value of 53 
is positive or negative (note that 53 cannot change its sign). The special feature of these two 
solutions is that in view of Eqs. ( |5o| ) and (^) three out of the six dominant susceptibilities are 
diverging with the same critical exponent. Namely 

XSDw ~ Xdsc ~ XCF ~ (A - Ac)"^ if 53 > 0, 

XCDw ^ Xssc ^ XSF ~(A-Ac)~^ if 53 < 0, (55) 
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where 7 = (\/5 — 2)/3 « 0.08. The divergence of the susceptibihties is thus weak compared to 
the mean field behavior x ~ (2^ ~ Tc)~^. 

However, for some initial conditions the solutions of Eq. (^) are not diverging but fiow 
towards the trivial fixed point 51 = (?2 = <73 = 54 = 0. This means that the RG flow of the 
couplings between the two saddle point patches does not develop an instability. In this case 
the restriction to the saddle-point patches is clearly insufficient and the low energy behavior is 
controlled by the remaining part of the Fermi surface. 



6.2 Towards functional renormalization 

The hypothesis of constant couplings gi, . . . ,g4 neglects the fact that in reality all these pa- 
rameters are functions of incoming and outgoing momenta ki, . . . , k4, three of which can move 
freely on the Fermi surface within the saddle point patches. Unfortunately, a true functional 
renormalization is presently beyond the reach of an analytical approach. Therefore we mimic 
the true momentum dependence by introducing, in addition to the constants gi, . . . ,g4 rep- 
resenting general values of the momenta (in the patches), other couplings corresponding to 
specific combinations of momenta. Thus we allow 53 (ki, . . . ,k4) to take four different values: 

ff3^^^ if kl + k2 = 

n 1 ^ J 93 if k3 - ki = Q 

33(ki,...,k4)^ <; if k3-k2 = Q 



53 



if |k3-k2-Q|, |k3-ki-Q|, |ki +k2| > O(v^) 



Similarly, gi (ki , . . . , k4) = takes the values gf or gi , 32 (ki , . . . , k4) = gf or 32 and (/4(ki , . . . , k4) 
gf'^^ or 514. We thus separate the couplings g^^^ ,gi'~'^ with zero total momentum, gf, g^ 
with a direct momentum transfer equal to Q and > 93 with an exchanged momentum transfer 
of Q from the general ones {gi, . . . , (^4), where none of these special relations among the in- and 
outgoing momenta applie^ In the transition domain (where for example < |ki -|-k2| < VA) 
the function is unknown, but, as we will argue shortly, its knowledge is not essential. 

This parameterization of the coupling functions is a natural choice given the structure of 
the one-loop RG equations (^9|) because the special couplings g^^^ ,g^ and g'^ get the strongest 
contributions from the PP, PHI and PH2 diagrams, respectively. We say that g^'^^ is resonant 
in the PP channel, g^ in the PHI channel and g'' in the PH2 channel. Furthermore the diagram 
with the largest contribution, say PP for a BCS coupling, again only includes BCS couplings 
and does not mix with non-BCS processes. 

For example, let k £ Pi and k' £ P2 such that (;(k, — k, — k', k') ~ ffif^^ and consider the 
RG equation for this process: 

-^gf = PP + PHI + PH2 (57) 

The dominant diagram PP, shown in Fig. |^, involves the couplings g^'^^ {k., p) and g^'^^ip, k'). 
Within our approximation they are replaced by the constants gz'~^^ or p^*^"^, respectively. The 
contribution from diagram PP to Eq. (^^ becomes 

PP = 2<p(A, 0) gi^^f''^, (58) 

where i3^(A, k) is the p-p bubble restricted to a saddle point patch of size p, given by Eq. ( p^ ) 
for t' = 0. 



^The special couplings are related to the couplings introduced in Section ^ by g^^ = 2gf — gf^ 
gBCS = gBCS^ = - g^ gi, = -g^ ggOS = ^SOS. 
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k' 

PP: p , 
k 




Figure 7: Two diagrams involved in Eq. (pT\). k and k' are typical vectors belonging to the patch 
Pi and P2, respectively; p is the integration variable 



Q(A,k + k ) := (60) 



By contrast, the non-resonant diagram PHI (also shown in Fig. ^ involves couplings like 
p(k, p — k — k',— k',p) which for almost every value of p do not satisfy one of the special 
relations ki + k2 = 0(\/A), ks - k2 = Q + 0(\/A) or ks - ki = Q + 0(\/A). This diagram 
therefore includes only the general couplings (72 and g-^ and no special couplings g^^^,g'' or 
g^. Thus 

PHI = 2Bp^,(A,k + k') 5253, (59) 
where i3^(A, k) — —Bpp{A, k — Q). It follows that the RG flow depends on the ratio _B^^(A, k-|- 
k')/_Bj^(A, 0). In order to obtain a closed set of equations we replace this ratio by a constant. 
Let us first define 

Bg, (A,k + k'- Q) 

which varies in principle between and 1. Its value is 1 if k + k' = Q, i.e. if the process to 
be renormalized is at the same time a g^^^ and a g^ . Such processes exist of course, but they 
will not infiuence the RG equations in a relevant way. According to Eq. ( p9| ) a region in the 
Brillouin zone of width ~ \/A or smaller can be safely ignored within logarithmic precision. 
We thus assume that k + k' — Q \/A or bigger. The biggest values are obtained if k + k' is 
parallel to the Fermi surface k + k' — Q = (k, k). For k > \/A we get from Eq. ( ^ ) 

q(A, k + k ) < — > -. (61) 

log^ ^ 

In the following we replace Q(A,k + k') by a constant a < 1/2. All the nearly resonant 
diagrams are treated in the same way, i.e. they include general couplings gi, . . . g4 only and 
their amplitude is reduced with respect to the resonant ones by a factor a < 1/2. 

Our approximation scheme leads to a set of RG equations for the special couplings gf, g2 , 93^'^ , si, 93 , 9?' 

d d ^ d/ X d\ , r\ d/ X d\ 

^Si = 231(32 -5i) + 2^3(53 - 53 ) 
^52 = (92) +(53) 

d BCS n BCS BCS , o /n ^ 

^53 = -2ff3 9a + 2aff3(2g2 - ffi) 

^53 = 2^2 53 + 2a 53 (32 -31 -34) (62) 
^S3 = '2{-2gigi + gzgt + gigs) + 2a 53(52 - 54) 

d BCS f BCS\'2 f BCS\2 

^54 = -(<73 ) -(<74 ) . 
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The general couplings gi , . . . , 54 are resonant in none of the three channels. The RG flow 
of (?i , . . . , 34 is thus given by Eqs. (^3|), with the right hand side multiplied by a. 

Eqs. (j^j can be rewritten in terms of the couplings which are associated with the dominant 
instabilities gi = -gi T 93, g± = - ± {2gi - g^,) and 5!^^ = gf""' ± gj'^^ 

= -{g±f ±'2aga{gi- g2+g4,) 

9± = -(Si)^ ±2033(31 +32 -34) (63) 



dl 

d BCS I BCS\'^ I r, /r, 

= -[9± ) ±2033(232-51 



We see that for a = it is a set of six independent equations, one for each instability. In 
fact, if the non-resonant diagrams are comple tely neglected, the RG becomes equivalent to 
the summation of ladder diagrams (see section S.4). Even for < a < 1 the special couplings 
associated with the different instabilities still do not influence each other, but each RG equation 
has a source term coming from the general couplings 31 , . . . , 34 . For a — 1 and initial conditions 
9i = 91, 92 = 92, gi = 93 = Sf^^ = gs, gf^^ = 34, Eq. (H) is recovered (since these 
conditions are then conserved by the RG flow). 

One can search for asymptotic solutions of the form g{l) = g-{lc — l)~^ of Eq. ( |63|) by solving 
the resulting algebraic equations for the 3. We first consider the possibility of diverging general 
couplings 31, ... ,34. In this case it follows from our analysis of Eq. (js^) that the asymptotic 
behavior of the general couplings is given by 31 = 0, 32 = —34 = l/{Goi) and 33 — ±\/5/(6a), 
depending on the sign of 33. By inserting this behavior into Eq. (p3) it is easily seen that a 
real solution for 3± , g% and 3±'^'^ requires a > ^80/81 ~ 0.994. But as we argued above, the 
appropriate values of a are < 1/2. 

It follows that for acceptable values of a a special coupling can only diverge if 31 = ... = 
34 = 0. This means that this special coupling constant diverges at a higher energy scale than 
the general couplings. 

The most striking difference to the "one-dimensional" solution is that one of the six couplings 
g%- can diverge, while all the others remain finite. This occurs here because the mixing of the 
flow for these couplings has been neglected on the basis of a phase space argument (see the 
discussion before Eq. (|5^)). This argument is certainly valid as long as the coupling function 
is slowly varying, but it may be questioned close to the instability, where the coupling function 
gets peaked. It is argued however in Section 6.4 that the non leading couplings can nevertheless 
stay finite in the case of an instability. 

The divergence of the leading coupling is characterized by g%. — —1. Eq. ( ^ ) then implies 
the asymptotic behavior x ~ (A — Ac)~^, corresponding to a mean field exponent. In view 
of Eq. (HTI) the coupling function 3*(k,k') diverges as (A — Ac)~^ everywhere on the Fermi 
surface and remains a smooth function upon renormalization, as anticipated in Section ^. 



6.3 Phase diagram 

For a cir cular Fermi surface the dominant instability is superconductivity, as discussed in Sec- 
tion [4.l| . For the square Fermi surface the flow equations ( |6^ ) show that there are several 
possible instabilities, s- and d-wave superconductivity, charge and spin density waves, charge 
and spin flux phases. The values of the initial couplings will determine which of these instabil- 
ities, if any, occurs flrst, i.e. at the largest energy scale. To be speciflc we consider an initial 
interaction consisting of on-site and nearest-neighbor terms 

Hi = U'^ n^rfiir ± X! ^ X! ' 
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Figure 8: The phase diagram for U > and a = 1/2. In the dotted region all the couplings are 
flowing to zero. The size of the dotted region depends on the parameter a, the other features are 
a-independent 



where Wr = Wor and Sr are the usual charge and spin operators on the lattice site r, the 
sum is over nearest neighbor bonds and U,V and J are parameters. The functional 

integral formulation of this model is of the form with a coupling function 

(?(fci, ...,k4)^U-{V^ J/4)ek3-k, + J/2 ek3-ki (65) 

where ek = — 2(cos kx + cos ky). 

The correct initial conditions for the flow equations would be the effective values of the 
couplings gi,. . . ,g4 at a cutoff Ao p^, when the flow enters the asymptotic regime. It is at 
present impossible to connect in a controlled way these effective low energy couplings to the 
parameters of the microscopic interaction because this would imply solving Eq. (^) without 
further expansion in powers of A. We therefore take the bare values of the couplings 

gi^gt^U-W-J, g2^g^ = U + 4V + J, 

93=91 = 93^93^''' ^U-W + 3J, g4 = g?'^'^ = U + 4V ~ 3J (66) 

as starting values instead of the (unknown) renormalized ones. This choice is a good approx- 
imation if U, V and J are small enough so that the couplings vary little before entering the 
asymptotic regime. 

We have solved Eq. 0631) numerically for these initial conditions and obtained the phase 
diagrams of Figs. |^ and ^for, respectively, positive and negative values of U. The predicted 
phase for the repulsive Hubbard model {U > V = J = 0) is a, spin density wave (SDW) as 
expected and quite strong nearest neighbor terms are needed to establish a flux phase or (only 
for attractive V) a d-wave superconductor. An unexpected feature of Fig. ^ is, that the SDW 
can be destabilized by positive values of J. 

In contrast to the "old" RG equations (|3|) our more elaborate scheme Eq. (^) produces 
only one diverging susceptibility while the others remain finite. Only at the phase boundaries, 
where the two neighboring phases are degenerate, both susceptibilities diverge. 

In a certain parameter range (the dotted region in Figs ^ and ^ all the couplings are flowing 
to zero. In this case the behavior is not necessarily dominated by the saddle points and the 
RG flow has to be followed to order A~^, including self energy corrections. 
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Figure 9: The phase diagram for U < and a — 1/2. 



6.4 A consistency test 



The approximation scheme presented in section 3.2 reduces to a generahzed random phase 
approximation (RPA) for a — 0. This approximation consists of summing a certain class of 
diagrams shown in Fig. px| . 

We have seen that although Eq. (^) at a > goes beyond the RPA, its solutions are 
asymptotically the same as for a — 0. The reason is that in our approximate treatment of the 
non-resonant diagrams only the general couplings gi , . . . , (74 intervene and not the (diverging) 
special couplings gf, §2 etc. It is not a priori clear whether this approximation is justified close 
to Ac. 

For example the (non-resonant) contribution of the diagram PHI of Fig. |^ to the renor- 
malization of g^^^ involves an integration over the coupling function (;(k, p — k — k', — k', p), 
where p is the integration variable moving along the one dimensional energy shell |^p| = A. It 
is an integral over a one dimensional curve in the space of momenta (ki,k2,ka). For special 
values of p, the curve crosses one of the planes specified by Eq. (^). For example for p = k' 
the coupling is equal to Qs^^ ■ The question is whether one can neglect the contribution close 
to this point even when g^'^'^ diverges. 

In order to answer this question we have to know the coupling function for small but finite 
total momentum. Within the generalized RPA one obtains 



BCS/i \ 

g± (A, q) = 



9± (Ao) 



l+pf^«(Ao)/;°Brp(A,q)dA 



(67) 



where q is the (small) total momentum|^and (;j?'^'^(Ao) is the (q- independent) initial value. Eq. 
( |67| ) can be obtained either by solviiig the RG equation ( |6^ for a = or by explicitly summing 

and 



the ladder diagrams shown in Fig. Ma). A similar expression is obtained for g'^ = ~g~ 
g'^ = 2g'^ — g^ as functions of the deviation of the momentum transfer from Q = (tt, tt). 
For negative values of g±^^ {Ko) the effective coupling diverges at a critical scale Ac. 
(A, q) has a maximum at q = 0, which diverges for A ^ Ac. 



For 



A > Ac, 



= p — k' in the example of section 3.2 
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Figure 10: The class of diagrams summed in the generalized RPA: a) particle-particle ladders and 
b) the particle-hole ladders. 



We have estimated the contribution of the peak in g±.'^^{A, q) to the non-resonant diagram 
by integrating this function over a curve in q-space. We find that such an integral diverges 
at most lilce log (A — Ac) and therefore is negligible compared to the resonant diagrams ~ 



(A — Ac) ^. We conclude that the approximation presented in section 6.2, namely evaluating 



the non-resonant diagrams with the constant general values gi, . . . , 34, is consistent. 

The weak divergence of the non resonant diagrams are also consistent with the fact stated in 
Section [s!^ , that non leading couplings and susceptibilities remain finite at Ac, since " dA log(A- 
Ac) < 00. Note that this behavior can only be obtained because the momenta are allowed to 
move continuously on the Fermi surface. If instead we would discretize the Fermi surface and 
replace the continuous coupling function by a finite set of constants, a divergence of one cou- 
pling Qc at a scale Ac would imply the divergence, at the same scale Ac, of all these couplings 
that have Qc appearing in the right hand side of the RG equation. 

Recently, the discretized RG has been studied in detail for a simpler model without van 
Hove singularities and without Umklapp scattering It was found that there is a factor 
between the biggest non-dominant couplings and the dominant ones, where A'^ is the number of 
patches. This is consistent with our result in the continuous case (A'^ 00), that non-dominant 
couplings stay finite at Ac while the dominant couplings diverge. Similar results were also found 
in the large- A'^ limit of half-filled A^'-leg ladders 



7 Special symmetries 

We will now discuss the phase diagram at half filling (Fig. ^ and P) in terms of special 
symmetries which turn out to be present on the lines separating two diff'erent phases. For that 
purpose it is convenient to write the model in the Hamiltonian formalism 

H = Ho + Hi- H' (68) 

with 

Ho — ^ ^ ^kflcrk, 
CT,k 
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Hi - 2l7 X] ^ki+k2,k3-k4 

ki ■ • ■k4 



X5r(ki, ka, kg, ki) 2^ c^^^c^ ,i^^c^,k3Cak4 



and 




,P)) 



■n-o-k. 



The Hartree-Fock term //' has been included in order to keep H particle-hole symmetric (i.e. 
invariant with respect to the canonical transformation Ccrk — > c^,^_,_q). 

The symmetry group of the noninteracting Hamiltonian Ho is extremely large. From any 
function dk we can build an operator Nd — 1/2 dk rio-k that commutes with Ho and 
thus generates a continuous group of symmetry transformations exp{iaNd) ■ We find that 
commutes with the complete Hamiltonian H if and only if 

(dki + dka - dk3 - dk4) 5(ki, . . . ,k4) = Vki,...k4. (69) 

If d is a suitably chosen d-wave function, such a symmetry can relate s- and d-wave su- 
perconducting order parameters by [Nd,OsSc] = —Odsc and [Nd,Odsc] = —OsSC- In this 
situation on can transform OaSC into Odsc and vice versa by a symmetry operation. As a 
consequence the susceptibilities for s- and d-wave superconductivity must be exactly equal, 
provided condition (p9|) holds. The symmetry Nd relates in a similar way spin- or charge 
density waves to the corresponding flux phases. It might therefore control the transition lines 
SDW/spin flux, CDW/charge flux and sSC/dSC. 

Similarly the operators Sd = 1/2 ^, ^ dk clj^r^rcr'Ccr'u. relate the spin density wave to the 

charge flux phase and the charge density wave to the spin flux phase. They commute with H 
if the following two conditions hold 

(dki - dk2 + dk, - dk4) sr(ki, . . . ,k4) = Vki,...k4 , . 

(dki -Hdk2 -dk3 -dk4)(l-X)5(ki,...,k4) =0 Vki,...k4. ^ ' 

Another symmetry introduced by Lieb and then further investigated by Yang and 
Zhang is generated by the pseudo spin operator r/s = Sk C|Q_kC|k, where the function 
Sk satisfies SQ_k = Sk- It turns a s-wave superconductor into a charge density wave and a 
d-wave superconductor into a charge flux phase (and vice versa), rjs commutes with Ho because 
of the exact nesting (Cq-k + ^k = 0) and it commutes with the full Hamiltonian provided 

EpSpff(p,P,k,k) = SkEp(2-^)ff(k,P,P,k) Vk 
Ski(l - X)g{]ii, . . . ,k4) - Sk33(k3,k2,k4, ki) + Sk43(k4, ka, ks, ki) = 

V ki, . . . k4, where k := Q — k. 

Finally Zhang considered the operators fid = 1/2 ^, ^ dk Co-q-k (t^t")^^, Cc,/k con- 
necting a spin density wave to a d-wave superconductor and a spin flux phase to a s-wave 
superconductor. The symmetry condition is of the same form as (M) but with Skj replaced by 
a function dki that satisfies dq_k = — dk. 

It is in general difficult to satisfy the conditions ( |69| ) to ( [n| ) . For example they do not hold 
for the U ~ V ~ J interaction. The only exception is the pseudo spin symmetry rjs which is 
exact for V — and Sk = 1. However the restriction of the model to the two saddle point 
patches has more chance of being symmetric. We take Sk = 1 everywhere whereas d(k) — 1 
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for k G Pi and d(k) = — 1 for k G P2- For this simple choice the symmetry generators Nd, Sd, 
rja and 11^ together with the total spin- and charge operators form a so(6) © so(2) Lie algebra. 
The commutation relations of the symmetry generators and the relevant order parameters are 
listed in reference psf . 

We further assume that the (initial) coupling function g{ki, . . . , k4) takes only four different 
values gi, . . .g4 as in Eq. (|5^). The symmetry conditions are then gg — for A'^^, ffi = for 
Sd, <?2 + 34 = 2(/i for rjs and (/2 + 34 = for Hd- These hyper-planes in our four dimensional 
coupling space define exactly the transition planes of the phase diagram (shown in Figs. ^ and 
^ for gi, . . . ,g4 parametrized by U, V and J) . 

We have thus shown that the transition planes of the phase diagram are fixed by exact 
symmetries of the gi, . . . , 34- model. This is a strong indication that the phase diagram shown 
in Figs. ^ and ^ is the correct one at sufficiently weak coupling. Such a determination of an 
exact phase diagram by simple symmetry considerations was also possible for a one-dimensional 
system Q, 



In the more general situation of Section 6.2 with ten instead of four coupling constants the 
symmetry conditions read 

gi ^ gi ^ gi'^^ = ga for Nd, 

91 = ffi = ff3 - 33 = for bd, 

n c\ d X BC S r\ d X BC S r\ c /rroN 

231-32-54 = 2^1-32-54 =253-33-33 =0 for r)s, (72) 

X , BCS X BCS n c rr 

32 + 34 = 32 + 34 =33-33 =0 for Yld. 

These symmetries are respected by our approximate RG equation (^2|) for every value of the 
parameter a. In fact it is easy to show that if one of these conditions is satisfied at the initial 
scale Iq it remains to be so at any scale I. This is not completely trivial since an arbitrary 
approximation scheme might violate the symmetries of the model. 



8 Conclusion 

In summary, we have analyzed systematically the instabilities of weakly interacting electrons 
with a square Fermi surface. Besides s- and d-wave superconductivity we have identified com- 
mensurate density waves and flux phases both in the charge and spin sector as the dominant 
instabilities. The transition lines of the phase diagram are fixed by exact symmetries and 
therefore robust for various approximation schemes. 

On a technical level, we have found that the dominant RG flow of the effective interaction 
in the limit of small energies is controlled by scattering processes with momenta close to van 
Hove points. Nevertheless, the low energy effective action contains relevant couplings between 
electrons everywhere near the Fermi surface. The couplings far away from the saddle points 
diverge at the same critical energy and with the same power ~ (A — Ac)~^ as the couplings at 
the saddle points. From this point of view, there is no sign of a scenario with strong effective 
couplings at the saddle points and weak couplings on the remaining Fermi surface. 

We have shown that the RG equations, although strongly coupled at the initial stage, 
become decoupled in the asymptotic limit of small energies. Thus the asymptotic result turns 
out to be similar to that of a generalized random phase approximation (RPA). The decoupling 
arises because the effective coupling function is strongly enhanced only for special configurations 
of the external momenta, i.e. in a small region of k-space. Thus the RG flow generates long 
ranged interactions in real space. 

By contrast, in one dimension superconducting and density wave instabilities remain cou- 
pled. The low energy excitations of a one-dimensional electron gas are constrained to two 
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privileged points in k-space: the Fermi points. This special geometry allows for a strong 
mixing between the various interaction channels. The saddle points of the two-dimensional dis- 
persion have a different status. They are privileged only due to the diverging density of states, 
but the low energy excitations exist on the whole Fermi surface. Therefore typical external 
momentum configurations of an effective interaction are resonant in at most one channel. 

Our analysis has further shown that a discretization of the Fermi surface in terms of a 
finite number of patches can enhance artificially the coupling between the different scattering 
channels in the low energy regime. In fact, in the immerical studies |2^, ^ different 
susceptibilities are found to diverge at a single energy scale. By contrast, our results show that 
only the susceptibility of the dominant instability diverges. This behavior has been referred to 
as the "moving pole solution" by the Russian school | pO| . 

The decoupling of the RG equations admittedly has only been established to leading loga- 
rithmic order in the energy cutoff A. This is justified for weak bare couplings for the square 
Fermi surface. Subleading contributions have to be taken into account if the Fermi surface is 
not nested at the van Hove filling or if the (bare) interaction is not small enough. 

In order to estimate how small the bare interaction must be, we recall that our approach 
requires a patch around the saddle point small compared to the size of the Brillouin zone 
(p <^ tt/S). On the other hand, the density of states of this patch has to be big compared to 
the remaining part of the Brillouin zone (in view of Eqs. (^) and (^^ this means log4/p^ <C 
log4p'^/A). The bare interaction qa^ ~ U must be small enough such that 



(The factor 2 appears because there are two saddle point patches in the Brillouin zone.) If for 
each of the two "<C" signs above, a factor of 10 is introduced, this amounts to [/ ~ 0.02 1. A 
less stringent factor of 3 for each of the two inequalities would correspond to 17 ~ 0.6 1. 
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